Optimized pulse sequences for suppressing unwanted transitions in quantum systems 
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We investigate the nature of the pulse sequence so that unwanted transitions in quantum systems 
can be inhibited optimally. For this purpose we show that the sequence of pulses proposed by 
Uhrig [Phys. Rev. Lett. 98, 100504 (2007)] in the context of inhibition of environmental dephasing 
effects is optimal. We derive exact results for inhibiting the transitions and confirm the results 
numerically. We posit a very significant improvement by usage of the Uhrig sequence over an 
equidistant sequence in decoupling a quantum system from unwanted transitions. The physics of 
inhibition is the destructive interference between transition amplitudes before and after each pulse. 
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I. INTRODUCTION 

The power of quantum information processing stems 
from the uniquely quantum phenomena of superposition 
and entanglement. Quantum bits (qubits) are two-level 
systems with states we can associate with the and 1 
of classical bits, but the valuable superposition quality is 
additional infomation coded in phase information called 
coherence Coupling of a qubit to it's environment 

entangles the two sytems, which effects a partial mea- 
surement on the qubit and introduces phase errors, a 
process known as decoherence [3] [I]. In NMR terminol- 
ogy, this corresponds to energy-conserving T2-type decay 
[5]. Qubits are also subject to non-energy conserving 
amplitude errors, or T^-type decay. Dynamical decou- 
pling seeks to suppress errors with appropriate pertur- 
bations. Application of 7r-pulses to a qubit-bath system 
- dubbed quantum "bang-bang" decoupling - effects a 
sort of time reversal in the system, which creates inter- 
ference and partial cancellation of environmental effects 
[4][6][7]- Further work showed that random or concate- 
nated pulse sequences could offer improvement over a 
sequence spaced equally in time [5] [5] [TO]. In this con- 
text Uhrig has shown that an optimal pulse sequence 
exists for decoupling a spin- 1/2 particle from a bath of 
harmonic oscillators [H][H]. The improvement of this 
sequence for preserving coherence has also been verified 
experimentally [T3] [H] . For a time interval T and an N- 
pulse sequence, the Uhrig Dynamical Decoupling (UDD) 
sequence suppresses decoherence up to order T N+1 . Yang 
then showed that UDD suppresses decoherence up to or- 
der T N+1 for a more general qubit-bath Hamiltonian [T5j . 
Indeed, Yang shows that the UDD is optimal for sup- 
pressing T\- or T 2 -type decay- corresponding to either 
phase or amplitude errors. The ideas of Uhrig have been 
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extended to protection of entanglement in two or more 
qubits [IS] . 

The present paper considers the problem of un- 
wanted transitions (decay and excitation) between two 
atomic levels due to a perturbing field. We exam- 
ine a scheme which suppresses environmental effects by 
quasi-instantaneous phase changes on a two-level system 
[T7][IB]- By applying a sequence of ultrashort 2-7T pulses 
to one state of the two-level system through an auxilliary 
level, this state acquires a phase shift of tt with respect to 
the other (See Fig.l). The phase shift leads to destructive 
interference between evolution terms in the time intervals 
between the pulses. Note that unlike the traditional spin 
problems where tt pulses are applied between the two 
levels, here we apply 2-7T pulses between one of the levels 
and an auxilliary level. This idea has been implemented 
in experiments in which fullerene qubits are decoupled 
from nuclear Rabi oscillations [19]. The question that 
we investigate in this paper is: what is the optimum se- 
quence of pulses to suppress unwanted transitions? We 
will show that the UDD sequence is optimal for a two 
level transition in an external field in the semi-classical 
limit. 

The plan of the paper is the following. In Sec. II we 
prove the optimality of the Uhrig sequence using first- 
order perturbation theory in the case of large detuning. 
In Sec. Ill we look at the exact case using the formalism 
of transfer matrices. In Sec. IV we prove again the opti- 
mality of the UDD for zero detuning. Sec. V reinforces 
these findings with numerical studies which suggest the 
optimality of the UDD across the all ranges of detunings. 
Finally, Sec. VI discusses applications and conclusions. 
In summary, we find that the UDD sequence offers a sig- 
nificant improvement over traditional equidistant pulse 
sequences for decoupling a two-level system from an un- 
wanted transition. 

Furthermore, we posit that using the UDD in experi- 
ments would improve the decoupling by a factor ~ ~ for 
resonant perturbations and by many orders of magnitude 
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FIG. 1. System under study: two coupled levels and an aux- 
iliary level 
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FIG. 2. Schematic diagram of pulse sequence. Notice that 
Si quantifies the timing of the i th pulse and r% is the time 
interval between pulses i and i — 1. That is, Tj = T(5j — 
Notice also that T = 2~Zil"i lr i> where i~i is a function of N 
such that T is fixed. 



two-level system and the field we seek to suppress. Of 
course, since the Schrodinger picture and interaction (or 
Dirac) picture are related by a unitary transformation, 
real observables or expectation values are unchanged. 
Hence, we begin with Eqn.[l]and write it as the sum of an 
unperturbed term, Hq = ^f-& z , and a perturbation, or in- 
teraction term, iJ 7 (i) = e -i<y *c7 + + h^e lu>t a-. Then 
the Hamiltonian becomes in the interaction picture: 

Hi(t) = e bHot ' h H I {t)e- lHot/h 

m 



e lAt +a-e~ lAt ] 



(2) 



where A = ui — ui' . 

After transforming the state vector into the interac- 
tion picture using the transformation I'F) ~ e lHot l h \%>) 
and rewriting the Schrodinger equation in matrix form 
in terms of expansion coefficients in the interaction pic- 
ture c e and c g , we find the coupled first order differential 
equations: 



in 



-g\-j 2 

in 

~2 



C~e(t) = -^ A %(t) 



(3) 



for increasingly detuned perturbations. For example, in 
Figure 3 we reproduce numerically some data from an 
experiment by Morton, where the equidistant sequence 
was used [19]. In addition, we have generated the data 
under the influence of the UDD sequence to illustrate 
graphically the improvement offered by such a sequence. 



II. LARGE DETUNING 

We examine a two level system, with levels labelled \g) 
and |e) coupled by a field which may or may not be on 
resonance, and real interaction parameter f2. We use a 
sequence of ideal 2ir pulses at times < 8\T < 82T < 
••• < 6{T < ••• < SnT < T coupling \g) to an auxilliary 
level \a) in order to effect a 7r-phase change in jg) (Fig- 
ures 1 and 2). We will see that because of the coupling 
between the levels, this phase change leads to destruc- 
tive interference between terms in the time evolution of 
the transition amplitude. We examine the time evolution 
of the transition amplitude for a general pulse sequence 
characterized by 5i, and ask for the optimum sequence. 

We use the semi-classical interaction Hamiltonian for 
a field of frequence uj' applied to a two-level system: 
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H(t) = —a z 
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Under the approximation of large detuning, A >> £7, 
we can use the machinery of first-order time-dependent 
perturbation theory. 



c e (t) - c e (t ) = Ac e (t,t Q ) 

n 



in_ . . 

-y c 9 (*o) 



lAt 'dt' 



2A 



( e iA(t-i ) _ 1)^0^) ( 4 ) 



For the remainder of this section, we set c g (0) = 1, 
which corresponds to an atomic population intially in 
the ground state. We can immediately reproduce the 
standard result for the probability of transition without 
a perturbing pulse sequence [3]: 



(5) 



Note that we can see the mechanical effect of the 7r 
phase shift on c g ; it sends Ac e (i,io) ~> — Ac e (£, £rj)- To 
introduce the effects of a pulse sequence, let = T(Sk — 
<5fc_i) be the k th time interval between the k th and k — l th 
pulses, and explicitly write out the evolution of the state 
vector. For times before the first 27r pulse is applied, the 
state vector evolves under the influence of the laser field 



as: 



t <t<to+n = |<7)+Ac e (£,i )|e) 



where a + = (a x + iu y ), o~ = [a x — i<Jy), and an 
interaction parameter, the Rabi frequency. 

In this section and Sec. IV, we will work in the inter- 
action picture, because it is the interaction between the 



The first pulse is applied after t\ has elapsed, effecting a 
7r phase change in the ground state, 

t = t +T! |*(t)) =-|ff)+Ac B (*o+Ti,t )|e) 
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after which the system evolves under a sort of time re- 
versal due to the phase change of \g), 

for t + Tl <t < t + Tl + T 2 
= -Iff) + Ac e (t +Ti,i )|e) - Ac e (t,t + n)\e) 

until the next pulse is applied. 

t =t + n + r 2 
|*(t)) =\g)+Ac e (t +r 1 ,t )\e) 

- Ac e (t + n + T 2 ,t + ri)\e) 

In general, after iV-pulses: 

N+l 

T=t +J2 T k 

fe=i 

|tt(T)) -(-1) W |.9) + 

AT p+1 p 

A£e(to + J2 Tfc ' *° + E r *) l g ) ( 6 ) 



p=0 



k=l 



k=l 



c e (i) 

Hence, upon substituting the explicit form of Ac e {t,to)'. 



JV 



£e(T) = ~2A E(- 1 ) P ( eiATp+1 - l)^ AE - lTfc (7) 



p=0 



Upon making the change of variables = T(5k — Sk-i), 
Eqn. [7] can be put into the form: 



| Ce (T)| 2 = ^ [1 + (-1) W +V AT + 2 £(-l)V AT ^] 

P =i 

(8) 

One can show that |c e (T)| 2 for an equidistant pulse 
sequence, Si = i/(N + 1), is: 

fl 2 AT AT 

\c e (T)\ 2 = — tan 2 — — - cos 2 — N even (9) 
A 2N + 2 2 

O 2 AT AT 

| Ce (T)| 2 = ^ tan 2 sin 2 ?± N odd (10) 

A 2N + 2 2 K ' 

This corresponds to one of the key results of Agarwal et. 
al. [IT]- One can see immediately that a large number 
of pulses suppresses the transition probability. We now 
look for an optimum sequence Si such that the transi- 
tion probability is minimized. The expression inside the 
brackets of Eqn. [H] is identical to Eqn. 10 in Uhrig's pa- 
per [IT], and the author uses the N degrees of freedom 
afforded by N pulses to set the first ^-derivatives with 
respect to T of Eqn. [8] equal to zero. The result is a sys- 
tem of coupled equations which can be solved for Si to 
yield: 



Si = sin 



2N + 2 



(11) 



Plugging this expression for Si into Eqn. gives an ana- 
lytic expression for the probability of transition: 



\Ce{T)? 



n 

2A 



E (-1: 



j e (iAT/2)cos(7rj7(n+l)) 



j=-n-l 



m 2 

~A 



2 (N + l) 2 J 2 N+1 (AT/2) 



(12) 
(13) 



The approximation follows from the Jacobi- Anger expan- 
sion, evaluation of a geometric series, evaluation of inde- 
terminate terms as a limit, and finally approximating a 
series by its leading term. Note that the approximation 
is good for AT < 2iV+2 [TT]. Figure 4 illustrates graphi- 
cally the efficacy of the UDD compared to the equidistant 
sequence for a relatively small pulse number. 



III. EXACT SOLUTION 

We now solve the Rabi problem exactly in the rotat- 
ing frame using the formalism of transfer matrices. Of 
course, since the Schrodinger picture and rotating frame 
are related by a unitary transformation, real observables 
or expectation values are unchanged. 

To begin, we use properties of Pauli matrices to rewrite 
the Rabi Hamiltonian, Eqn. [TJ 



H(t) 



h'jj 



cr,e 



(14) 



To eliminate the time dependence of the Hamiltonian, we 
switch to a rotating frame. 



W) 



(15) 



Under this transformation the Schrodinger equation in 
the rotating frame becomes: 

^ / ) = h Y°^') + h Y°^') ( 16 ) 

By writing our state vector in the two-level basis in the 
rotating frame, this becomes in matrix form: 



/ 1 -A il \ 1 <' 
Q A 



(17) 



Eqn. 17 is a system of coupled first order linear differ- 
ential equations which can be solved to yield a generalized 
Rabi frequency, fl R ; = \J A 2 + f2 2 and the time-evolution 
of the coefficients of expansion in the rotating frame. 



<(*) 

<£(*) 



M(t) 



4(0) 

4(0) 



(18) 



where M(t) is the unitary matrix: 



in 



sm(±fl R t) 



^sin(i^) 

A 



cos(ifi fl t) + fA S i n an R t)n 
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As a mathematical note, we can now see why the ro- 
tating frame was the best choice for solving the prob- 
lem exactly. We use the rotating frame because the SU 2 
group property has a desirable form when we transform 
to a picture in which the Hamiltonian is time indepen- 
dent; namely, for M(r), N(t) £ SU 2 , M{ti)N{t 2 ) = 
MN(ti +t 2 ). This is desirable because in physical terms 
this says that the evolution for time t\ followed by evo- 
lution for time r 2 is the same as evolution for total time 
Ti + t 2 . 

Now we define a matrix which represents the appli- 
cation of an ideal 27r-pulse to the ground state through 
the auxilliary level. Each of the two-levels is an eigen- 
state of the spin operator a z , so each state possesses the 
unique spinor property that a rotation by 2tt produces a 
ir phase shift [5]- Hence, a 27r-pulse changes the sign of 
the ground state. 



Again, we describe pulses by matrices, II, and we form 



an equivalent expression to Eqn. 20 in the interaction 
picture. 



Ce{T) 



= M(T N+1 )f[TlM(n) f^jjj 



(25) 



We can actually find a simple analytical solution to 
Eqn. [25] quite easily. The key is to simply starts multi- 
plying out some of these matrices and use the trigono- 
metric identities sin a cos /3±sin j3 cos a — sin (a ± f$) and 
cos a cos /3 =F sin a cos f3 — cos (a ± Upon making the 
notational definition 0,; = — 7j_i + • • ■ + (— l)Vi one 
can see by induction that Eqn. |25| becomes: 



Ce(T)\ _ ( {-l) N cos{\9,@ N+ i) -isva(\Q.Q N+li 



i{-l) N sm{\VLQ N+l ) cos{\ttQ N+l ) 



n 



1 

-1 



(19) 



Ce(0) 
Cg(P) 



(26) 



Using this formalism, we can easily write down an ex- 
pression for c' g and c' e after an ./V-pulse sequence. 



c'e(T) 



= M(r N+1 )UM{r N )...IlM{n) 



JY 



= M(t n+1 ) Y[UM(ri) 



c'M 

^(0) 



Therefore, solving finally for the excited and ground 
state probabilities, we have: 

\c g {T)\ 2 = |c 9 (0)| 2 cos 2 (iQe^+i) + |c e (0)| 2 sin 2 (i09 w+1 ) 

(27) 

|Ce(T)| 2 = |C 



| 2 cos 2 (^e A r +1 )-|-| C3 (0)| 2 sin 2 (ifie 



(20) 



N+lJ 

(28) 



In order to decouple the levels from Rabi oscillations, 



In Sec. V we turn to the numerical study of Eqn. [20] for we want to trap the state such that |e 9 (T)| 



different pulse sequences, and hence different expressions 
Ti. Numerical studies , such as Figure 5, suggest that the 
UDD is optimal over all ranges in detunings. 



IV. SMALL DETUNING 

We now examine the limit of zero-detuning, or equiv- 
alently, large pulse numbers.. If we say that r w T/N, 
this case is the limit ^? ps 0. We work in the interac- 
tion picture and solve the coupled differential equations 
in Eqn. [3] after setting the product At = and find, upon 
solving for coefficients in the interaction picture: 



and |c e (T)| 2 w |c e (0)| 2 . Hence, we seek a pulse sequence 
such that Qn+i is minimized. 

Our expression for Oat+i can be put into a more useful 
form by the change of variables n = T(Si — We 
find that: 



N 



e N+1 - t(i + 2(-i) N+1 5](-i) fe 4 



(29) 



fc=i 



To see how well the equidistant sequence does in decou- 

and use the arithmetic 



pling the levels, we set Si 
identities: 



N+l 



c g (t) = c g (to) cos(ifit) - ic e {t ) sin(iQt) (21) 
c e (t) = c e (to)cos(|Ot) - icg{t )sm{\Slt) (22) 

In the spirit of the previous section, we define another 
unitary transfer matrix M(t): 

M(t)=( COS ^ f ~ is %f) (23) 
\—ism(±ilt) cos(|iM) J 

Hence, the time evolution of c g and c e between to and 
t is give by: 



k=l 
m 

E 

k=l 



(-l) k k = 



(-l) k k 



rn 



(m+1) 



m even 



m odd. 



We find that: 



Ojv+i 

@JV+1 



T 



N 




N even 
N odd 



(30) 
(31) 



(32) 
(33) 



Ceit) 
Cg(t) 



= M(t) 



c e (to) 

Cg{t ) 



(24) 



Hence, we can decouple the levels to arbitrary precision 
given a large number of pulses. In fact, any number of 
odd pulses will decouple the two levels completely in the 
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sense that \c g {T)\ 2 = \c g (0)\ 2 and |c e (T)| 5 



|c ff (0)| : 



This is easy to understand in terms of the physical effects 
of the 7r phase shift on the ground state. The pulses 
effects a sort of time reversal where the system retraces 
it's steps after the pulse is applied. Hence, after an odd 
number of pulses, the system ends up exactly where it 
started. Then why is the Uhrig pulse sequence optimum 
for a general qubit-bath system? This is because the time 
reversal effected on the system affects only the qubit. 
Hence, the qubit tries to retrace it's steps, but the bath 
is still constantly reacting to the instantaneous state of 
the qubit. Since most of the time we assume that the 
initial qubit-bath state is a product state, after any finite 
time interval the system irreversibly becomes a mixed 
state. Therefore, the time reversal is imperfect, and the 
equidistant sequence is not optimum. Then the question 
arises: what is the optimum sequence for any TV in our 
bathless system? 

To answer this question, we seek a function Si such 
that: 



\N+1 



N 

k=l 



(34) 



Which is simply setting Eqn.[29| equal to zero, dividing 
out nonzero T, and multiplying both sides by (— l) w+1 . 
Upon inspection, one can see that the left hand side is 
actually the first derivative of Eqn. [8] at AT = 0. By 
construction, the UDD satisfies the hP* derivative equal 
to zero. Specifically, the UDD satisfies the first derivative 
equal to zero. Hence, the UDD again optimizes the de- 
coupling of the levels in the sense that the system returns 
to it's initial state. 
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FIG. 3. Plot of the evolution of the ground state probability 
\c g \ 2 versus time for 8 pulses and A/fl = 0. Note that the 
Uhrig sequence is more effective than the equidistant sequence 
in preserving the intial state of the sytem at t = 1.5 since 
| Cg | 2 = 1 only for the Uhrig sequence. 



V. NUMERICAL RESULTS 

We now examine the efficacy of different sequences for 
preserving the initial state of the system, that is, how 
close |c e (T)| 2 fa |c e (0)| 2 . The pulse sequences to be stud- 
ied are: (i) no pulses, (ii) equidistant sequence, Si = ^Vj 

and (iii) UDD sequence, Si — sin 2 ( 2 £ l +2 ). The following 
plots illustrate the efficacy of each sequence for suppress- 
ing a transition to the excited state when the popula- 
tion is initially in the ground state. Figure 3 shows Rabi 
oscillations as a function of time, and how they are af- 
fected by the pulse sequences. The pulses effect a per- 
fect time reversal on the system for zero detuning, and 
the UDD has a symmetry which always returns the sys- 
tem to where it started. In contrast, the efficacy of the 
equidistant sequence depends on the parity of the pulse 
number. Note that both sequences have some symmetry, 
which the equidistant sequence loses with non-zero de- 
tuning, causing the UDD to be more effective for even or 
odd pulse numbers. Compare this figure to experimental 
data by Morton [T^]. Figure 4 compares the evolution 
of the system for different pulse sequences for a large 
detuning and small pulse number. Here the equidistant 
sequence offers a suppression of the transition by a fac- 
tor of about 10 2 , and the UDD sequence suppresses the 
transition by a factor of 10 5 . Note how the UDD causes 
the transition probability to remain relatively flat for a 
period of time; this is the result of the first ./V-derivatives 
being set to zero. 
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FIG. 4. For n = 5 and A = 10S7, pictured is a comparison 
of transitions probabilities for no pulse sequence, equidistant 
sequence, and UDD sequence, respectively, as a function of 
time. 

Finally, Fig. 5 compares the efficacy of the sequences 
over a range of pulse numbers, plotted with the ana- 
lytic solutions for the large detuning case we derived. 
It seems that the analytic solution for large detuning is 
still a good fit for this intermediate detuning. The Uhrig 
sequence is more effective than the equidistant for sup- 
pressing the transition for all n > 1. For n = 1, the 
UDD and equidistant sequences give the same result be- 
cause the UDD degenerates to the equidistant sequence 
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FIG. 5. Plot of the probability of transition |c e | 2 versus pulse 
number for A/fl = 1 for 2-11 pulses. 



for n — 1. Note that each of the pulse sequences have 
been multiplied by an arbitrary factor of 20 to make the 
graph more readable. The transition probability with 
the UDD sequence drops very rapidly, reaching ~ 10 , 
- lCT'V 10~ 10 , and - 1CT 14 , for n = 2, 4, 6, and 8, 
respectively, compared to ~ 10~ 2 for the equidistant se- 
quence across this range of pulse numbers. As a note, 
in Fig. 5 the system is allowed to evolve through half a 
Rabi cycle, that is T = 7r/0^, so the transition proba- 
bility under no perturbing influence is a maximum for a 
given A/O. In summary, the Uhrig pulse sequence proves 
to be many of orders of magnitude more effective for all 
values of A/O. 



VI. CONCLUSION 

The decoupling of the levels of a two-level system has 
recently been realized in the experiments of Morton [TDJ . 
In this experiment, nuclear fullerene qubits are decou- 
pled from nuclear Rabi oscillations by taking advantage 
of the coupling between a nuclear and electronic qubit. 
Through a repeated application of a very fast phase gate 
a nuclear spin qubit was bang-bang decoupled from a 
permanent driving field. The speed of the phase gate 
and the resonance condition means that the equidistant 
sequence works very well in decoupling the system. How- 
ever, as we have shown, the UDD sequence offers an im- 
provement over the equidistant sequence for any value of 
the detuning. Hence, this decoupling scheme works well 
in tailoring the interactions in quantum systems, lend- 
ing itself to broad applicability in quantum computing 
schemes [2D]. In addition to decoupling a system from 
its environment, the effect is required in some quantum 
computing schemes [21] . 

In conclusion, we have proved analytically that the 
Uhrig Dynamical Decoupling pulse sequence optimizes 
the inhibition of unwanted transitions in a two-level sys- 
tem for the case of large detuning as well as small de- 
tuning. We have further confirmed our analytical results 
numerically. These results may prove fruitful for con- 
struction of a quantum computer, where suppression of 
decoherence and unwanted transitions arc important con- 
siderations. 
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